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ABSTRACT 

We present deep spectroscopic measurements of 18 distant field galaxies identified 
as gravitationally-lensed arcs in a Hubble Space Telescope image of the cluster Abell 
2218. Redshifts of these objects were predicted by Kneib et al. (1996) using a lens- 
ing analysis constrained by the properties of two bright arcs of known redshift and 
other multiply- imaged sources. The new spectroscopic identifications were obtained 
using long exposures with the LDSS-2 spectrograph on the William Herschel Tele- 
scope and demonstrate the capability of that instrument to new limits, R ~24; the 
lensing magnification implies true source magnitudes as faint as R ~25. Statistically, 
our measured redshifts are in excellent agreement with those predicted from Kneib 
et al.'s lensing analysis which gives considerable support to the redshift distribution 
derived by the lensing inversion method for the more numerous and fainter arclets 
extending to R ~ 25.5. We explore the remaining uncertainties arising from both 
the mass distribution in the central regions of Abell 2218 and the inversion method 
itself, and conclude that the mean redshift of the faint field population at R ~25.5 
(B ~ 26-27) is low, <z>=0.8-l. We discuss this result in the context of redshift dis- 
tributions estimated from multi-colour photometry. Although such comparisons are 
not straightforward, we suggest that photometric techniques may achieve a reason- 
able level of agreement particularly when they include near-infrared photometry with 
discriminatory capabilities in the 1 < z <2 range. 
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1 INTRODUCTION 

Gravitational lensing is now established as a highly- 
successful technique for studying both the mass distribu- 
tion on various scales in rich clusters at intermediate depth, 
and for constraining the properties and distances of very 
faint sources magnified serendipitously via the lensing clus- 
ter. Smail et al. (1994, 1995a) attempted to break the degen- 
eracy between these two important applications by studying 
the lensing signal in a number of rich clusters chosen to lie 
at different redshifts. Although they concluded the bulk of a 
faint galaxy sample limited at I ~25 has a redshift distribu- 



* Based on observations obtained on the William Herschel Tele- 
scope at the Observatorio del Roque de los Muchachos, La Palma. 



tion with a median close to one, their conclusion depended 
somewhat on the mass distribution in the more distant and 
hence less well-understood clusters. 

The discovery of multiply-imaged systems in the cores 
of rich clusters (Mellier, Fort & Kneib 1993, Kneib et al. 
1993, Smail et al. 1995b) offers an alternative and more 
promising route. The mass models for clusters containing 
such multiple images are sufficiently constrained by the lo- 
cation and orientation of the images that it is possible to 
use these models to match the predicted shear of the fainter 
images with their observed shapes directly. In the case of 
multiple images, the requirement that each image must come 
from the same object in the source plane provides many con- 
straints leading to a rather precise redshift prediction. Singly 
imaged arclets generate fewer constraints but nonetheless, 
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a most likely redshift may be determined for each arclet, 
and useful statistical results can be obtained for the back- 
ground redshift distribution, provided the sample size is 
large enough. A first application of this so-called 'lens inver- 
sion' technique was discussed by Kneib et al. (1994) for the 
cluster Abell 370. Individual redshifts were obtained for 30 
candidate arclets limited at B—27 with axial ratios a/b >1.4 
measured on ground-based images. 

The inversion technique was considered in more de- 
tail by Kneib et al (1996, hereafter KESCS) on the basis 
of an impressive image of Abell 2218 (z=0.175) obtained 
with the refurbished Hubble Space Telescope (HST), Figure 
|l| The important role of HST imaging in lensing studies 
was demonstrated in two respects. Firstly, via the superla- 
tive resolution of HST, the multiply imaged nature of several 
of the most strongly lensed arcs was revealed in spectacu- 
lar detail. Secondly, in the absence of ground-based seeing, 
the fainter sheared images were individually identified with 
much greater confidence than was possible using ground- 
based data. This, in turn, enabled accurate inversion for a 
greater number of faint sources. Predicted redshifts were de- 
termined by KESCS for ~80 arclets in Abell 2218 to i?=25 
and a formalism was developed to predict the uncertainties 
arising in the inversion from a number of sources. The limi- 
tations of lensing inversion as a general surveying technique 
were also discussed. 

Although impressive progress has been made in locating 
high redshift galaxies via deep (K <20) magnitude-limited 
redshift surveys (Cowie et al. 1995, 1996) and via various 
colour selection criteria such as those sensitive to the pres- 
ence of a Lyman limit in the near ultraviolet (Steidel et 
al. 1996), an important advantage of lensing inversion as 
a purely geometric technique over other methods is that it 
does not require the presence of a particular spectral fea- 
ture within some restricted observational passband. This 
suggests it should complement the other techniques by pro- 
viding an accurate mean redshift distribution at various lim- 
its even if, ultimately, inversion through several clusters at 
various redshifts may be required to secure robust results. 

Given the promise of the method it is clearly important 
to investigate the uncertainties of lensing inversion by direct 
spectroscopic measurement of the brighter arclets. This pa- 
per is concerned with executing such a test for Abell 2218 
whose mass model is highly constrained from 7 multiple im- 
age sets, two of which already have spectroscopic redshifts 
(Pello et al. 1992). Whereas only a few of the brighter arclets 
are amenable to direct spectroscopy with current facilities, 
the test should still be valuable since inversion is primarily 
a geometric technique. The estimated redshift for a particu- 
lar arclet relies on the cluster mass distribution, the relative 
distances of the source and lens, and uncertainties in the 
source and image shape (whose contribution can be readily 
accounted for - see KESCS). Thus a comparison of the spec- 
troscopic and predicted redshifts for the brighter arclets can 
be made as valid a test of the method as that, more difficult, 
test applied to the entire sample of arclets. 

We have begun a survey of those arclets amenable for 
spectroscopic study with the 4.2m William Herschel Tele- 
scope (WHT). Our aim is to verify the inversion technique 
and, if necessary, correct the residuals and produce an im- 
proved mass model. This in turn will lead to better inversion 
for the more numerous and fainter arclets. As part of this 



programme, one arclet in the sample (#384 in the num- 
bering scheme of Le Borgne et al. 1992) has already been 
discussed by Ebbels et al. (1996) owing to its high redshift 
(2=2.515) and the spectacular agreement with the lensing 
prediction (z = 2. 8^0 3)1 making it an object of interest in 
its own right. A plan of the paper follows. In §2 we discuss 
the selection criteria for our sample of arclets in Abell 2218 
and present the journal of observations and spectroscopic 
results. In §3 we compare the observed redshifts with those 
predicted and consider the uncertainties. In §4 we review the 
mass model and illustrate those changes which could still be 
permitted whilst being consistent with the new sample of 
arclet redshifts. Finally, in §5 we return to the question of 
the mean redshift of the faint population following the dis- 
cussion of KESCS. §6 summarizes our conclusions. 



2 OBSERVATIONS AND SPECTROSCOPIC 
REDSHIFTS 

The primary source material for our arclet sample is the 
HST WFPC2 image of A2218 as published by KESCS and 
shown in Figure |l|. The three orbit exposure totaling 6500s 
was taken in the F702W passband and the magnitudes con- 
verted to standard R using corrections from Holtzman et al. 
(1995). The resulting object catalogue reaches R — 26 with 
a completeness of 55% at R ~ 25. Ground based (B — I) 
colours for our targets were measured from images obtained 
with the COSMIC imaging spectrograph on the 5-m Hale 
telescope at Palomar during June 1994 and June 1995. A 
thick 2048 2 TEK CCD (0.284 arcsec pix~ v ) was used and 
individual exposures of ~ 500-1000s were taken using in- 
field dithering (by ~ 60 arcsec). The images were reduced in 
a standard manner with IRAF using both twilight and on- 
source flat fields and including cosmic ray rejection. Stacked 
exposure times totaled 16.5 ks in B and 21.7 ks in 7, reaching 
la limiting surface brightnesses of \jlb ~ 28.2 and hi = 26.8 
mag arcsec -2 , while the final seeing on the frames measured 
1.20 arcsec in B and 0.95 arcsec in I. B & / magnitudes 
were measured from these deep exposures. Due to the faint- 
ness of these objects and the likelihood of contamination by 
nearby bright cluster galaxies, we measured these colours 
by hand. First, using IRAF's IMREPLACE task, a rect- 
angular region enclosing the target object was replaced by 
fitting a surface to the perimeter pixels and then using the 
distribution of the perimeter pixel values to add sky noise 
to this surface. Then a difference image was formed between 
the original and object-replaced images such that only the 
flux from each object remained. Magnitudes were measured 
within circular apertures of radius 1 or 2 arcseconds and 
colours obtained from measurements within the same aper- 
ture in each band. 

In selecting our arclet sample from the HST image, we 
employed three criteria. Firstly, we selected objects which 
looked distorted, indicating that they might be gravitational 
arclets. As in KESCS, we maximized the probability of se- 
lecting lensed images by restricting the sample to those im- 
ages whose orientation lies within 45 deg of the local shear 
direction as predicted by the mass model. Secondly, to max- 
imize the probability of securing a reliable spectrum with 
measurable features, we applied an integrated magnitude 
limit of B ~ 25 (corresponding to R ~ 24) and gave highest 
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priority to objects bluer than (B — I) ~ 2.2. The former con- 
straint corresponds to the faintest objects within reach of the 
LDSS-2 spectrograph on the 4.2m WHT and the colour cri- 
terion is about half a magnitude bluer than the cluster E/SO 
sequence and increases the probability of finding emission 
lines in the spectrum. Finally, we were mindful that spec- 
troscopic redshifts for some arclets would be more important 
in the verification of the inversion method. The uncertainty 
for each arclet is a function of location, shape and redshift 
and thus by selecting those arclets for which the expected 
error is small, we can provide a more stringent test of inver- 
sion. These criteria were applied to 572 objects detected at 
2 sigma above the sky in the WFPC2 image and produced 
a working catalogue of some 65 candidate arclets. 

All observations were conducted with the LDSS-2 multi- 
object spectrograph (Allington-Smith et al. 1994) on the 4.2 
metre William Herschel Telescope and are summarized in 
Table [j]. In designing multi-slit masks, we were constrained 
by the fact that most of the highly-magnified arcs lie close to 
the cluster core. This caused problems with crowding since 
we insisted on a minimum slit length of 10 arcseconds to 
enable good sky subtraction. A second obvious limitation is 
the orientation and curvature of the arclets. Although with 
ground-based seeing only the giant arcs are noticeably non- 
linear in form, the orientations of the larger arclets must be 
taken into account in order to maximize the amount of light 
captured by each slit. Each mask had one primary target of 
this type and thus we were able to base the mask orientation 
on the inclination of this large arclet alone. 

Accordingly, from the above sample, 8 multi-slit masks 
were constructed each containing between about 5 and 15 
arclet targets which allowed for some duplication between 
the two runs. In addition to the main targets, spectra were 
obtained for several other objects, for instance where a part 
of the slit happened to fall on a second object near the main 
one or where a suitable arclet candidate was not available. 
These 'serendipitous' objects are at present included with 
the main sample but nagged as such on the following dia- 
grams. In summary therefore, the spectroscopic sample con- 
sisted of 61 targets of which 37 were bona fide arclet can- 
didates. The colour magnitude diagram for our targets is 
presented in Figure ^, superimposed on that for the whole 
cluster field. Our magnitude and colour limits can clearly be 
seen as can the similar (B-I) colours of several of the multi- 
ple image systems. The observations were completed during 
May/June 1995 and June 1996 under mixed conditions. A 
total of 7 masks were exposed as shown in Table |l|. 

The data were reduced using the LEXT package (Col- 
less et al. 1990, Allington-Smith et al. 1994). Frames were 
debiassed and median combined taking into account small 
(~ 0.3 pixel hr" 1 ) shifts due to instrument flexure. Dis- 
persed tungsten flat fields were used to correct for along- 
slit variations in illumination and pixel to pixel sensitiv- 
ity changes with wavelength. Third order polynomial fits to 
CuAr arc lines provided wavelength calibration to a r.m.s. 
residual of ~ 2.0 ± 0.1 A. Sky subtraction was performed in 
two steps. Firstly we fit the summed slit profiles, exclud- 
ing the object rows and edge effects. The deviation from 
this overall profile was then fit column by column to remove 
wavelength variations in the profile. Finally the object rows 
were extracted with a Gaussian weighting to optimize sig- 
nal to noise (except for cases where the target was partially 
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Figure 2. The colour magnitude diagram for the arc sample su- 
perimposed on that for the whole cluster field. Open circles de- 
note non-sample objects which were not obviously cluster mem- 
bers and the dashed lines show our magnitude and colour criteria. 
Objects outside this region were selected primarily by morphol- 
ogy or are serendipitous. Note the similar (B — I) colours of the 
multiple image systems (384,468) and (359,328,389). 



merged with another object, e.g. #273). Redshifts were de- 
termined independently by two observers (RSE and TMDE) 
using Karl Glazebrook's REDSHIFT utility and are sum- 
marized in Table ^. We assigned a quality flag q to each 
spectrum according to the reliability of the redshift deter- 
mination. Redshifts determined from more than one strong 
line, or a single line and continuum with good signal to noise 
(S/N) were assigned q = 1, designated 'certain / probable'. 
Those determined from one or more noisy lines were assigned 
q = 2, designated 'possible', and those for which no redshift 
could be obtained were given q = 3. 

Of the 37 arclet candidates for which slits were cut, 7 
failed to yield spectra of adequate signal/noise (never reach- 
ing S/N con t — 3). There are several reasons for this. Firstly, 
a major limiting factor at these faint limits is the effective- 
ness of sky subtraction. This is largely influenced by the sur- 
face brightness of the target compared to that of the night 
sky (nominally fin, s kv — 20.4 mag arcsec~ 2 at La Palma in 
dark time.) Figure shows the success rate of our survey 
with both integrated R magnitude and surface brightness, 
Hr. It is evident that our sample is strictly limited at the 
/in = 24 level and that our percentage success rate for deter- 
mining firm redshifts decreases as the mean surface bright- 
ness falls (100%, 54% & 28% in the three bins from fj, R = 22 
to hr = 24 respectively.) 6 of the 7 low S/N arclets have 
surface brightnesses fj,R > 23.2 (< 7% of the night sky) and 
integrated magnitudes of B > 24.0 (R > 22.0), the faintest 
limit attained by previous LDSS-2 surveys. In addition, in 
all but one case, less than 8 arcseconds of sky was avail- 
able for sky subtraction due either to contamination by a 
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Figure 1. HST WFPC2 image of the cluster Abell 2218 in F702W showing the sample of objects selected for spectroscopy (listed in 
tabic ^). Those for which we have obtained good redshifts are marked with filled labels. 

Table 1. Spectroscopic Observations. The LDSS-2 medium blue grism was used for all masks except mask 1 for which the medium 
red grism was used. Both grisms give a resolution of ~ 12A. During the 1995 run a thinned Tektronix 1024 2 CCD (TEK-1) was the 
detector, giving 0.59 arcsec pix~ 1 !k5.3Apix— 1 while a new Loral 2048 2 CCD (LOR-1) became available for the 1996 run, with a scale 
of 0.37 arcsec pix~ 1 &c3.3Apix~ 1 . Final stacked on-source exposure times were as follows: mask 1 - 19.8ksec, mask 2 - 19.1ksec, mask 
4 - 10.8ksec, mask 730 - 32.4ksec, mask 444 - 45.5ksec, mask323A - lO.Oksec, mask323B - 7.2ksec. NOTES: (1) Identification numbers 
follow Le Borgne, 1992. (2) Duplicate masks. These contained the same objects except that the bright target #323 on mask 323A was 
replaced on mask 323B by fainter targets to increase the number of spectra obtained. (3) Numbers in parentheses denote exposures in 
poor conditions. 



# Mask 


Targets 1 


Dates 


Total Exposure (ksec) 3 


Seeing(arcsec) 


1 


138,117,231,321,456,355,308,236,408 


1995 May 27 


10.8 


1.0 






1995 May 29 


9.0 


1.0-1.5 


2 


129,171,229,328,384,468,492 


1995 May 28 


5.4 


1-2 






1995 May 30 


11.9 


1.0 


4 


238,273,384,467 


1995 June 1 


10.8 


0.7-1.0 


730 


730,158,167,179,230,254,257,342,464,476,461 


1996 June 18 


14.4 


1.1 






1996 June 19 


7.2 (3.6) 


1.0 




a 


1996 June 20 


7.2 


1.0 


444 


444,117,145,236,242,256,262,297,306,307,308,344,381,431 


1996 June 16 


7.5 (7.5) 


0.8-1.3 






1996 June 17 


15.0 


0.8-1.0 






1996 June 21 


18.0 


0.9-1.1 


323A 2 


323,132,159,167,190,200,205,223,288,295,317,321,389 


1996 June 15 


10.0 


1.0-1.5 


323B 2 


277,132,159,167,190,200,205,223,288,295,257,321,389 


1996 June 20 


7.2 


0.9 



second galaxy in the slit or simply the length of the slit. Fi- 
nally, it is possible that misalignment of the slits from their 
required positions could have led to low signal to noise in 
some spectra. To investigate this, we took all the targets 
with 22 < R < 23 and measured the average count level 
in each spectrum in a ~ 1000A wide region spanning the 
position of the R filter (6500A). If most slits were aligned 
well with the target objects, there should be a correlation 
between the count level in the spectra and R magnitude. We 
performed this test for 15 arclets and indeed found a linear 
correlation of spectral counts with io~ when three out- 
liers were discarded. At least one of the outliers could be 
explained by poor slit orientation, suggesting that misalign- 
ment of the slits with object positions was not a significant 
cause of light loss in the majority of cases. 

Several spectra show good continuum signal to noise, 
but no recognizable features and it is reasonable to ask 
what redshift ranges are implied by the non-detection of 
features in these spectra. Following Colless et al. (1990), we 
can use those spectra where [Oil] emission was detected to 
calculate our sensitivity to detection of emission features. 
Given a continuum S/N value one may obtain a limit on 
the minimum equivalent width of emission which could re- 
sult in a detection for a particular spectrum. As the con- 
tinuum S/N increases, we should become increasingly sensi- 
tive to lower equivalent widths of emission. Figure ^| shows 
a plot of the continuum S/N adjacent to [Oil] versus the 
observed equivalent width of [Oil] for the 15 arclets from 
which [Oil] emission was detected. For an unresolved line, 
continuum S/N (S/N cont ), line equivalent width (W\ i0 b B ) 
and line S/N (S/Nu ne ) can be related to each other, via 
simple assumptions, through the relation S/N co „tW\ : obs ~ 



v2n&ii ne S/Nii ne . (We assume the line profile is a Gaussian 
of sigma <7u ne ~ 0A25FWHM. For our ~ 12A resolution, 
this gives on ne ~ 5.1.) Plotting this relation on the dia- 
gram shows that our detection limit corresponds to a value 
of S/N line ~ 5. 

We now proceed to use this detection limit to put 
bounds on the possible redshifts of those (g = 3) arclets 
with good S /N but no visible features. Given some signal to 
noise ratio for a piece of continuum, one may calculate using 
the relation of Figure ^ the corresponding minimum equiva- 
lent width of emission which could have been detected there. 
Combining this limit with the distribution of [Oil] equiva- 
lent widths seen in field galaxies from other surveys, we may 
derive a probability that there could be [Oil] emission at 
this point without it having been detected. The presence of 
[Oil] at this point defines a redshift and so this translates 
into a probability Ps/n( z ) that the galaxy could have this 
particular redshift without any [Oil] being detected. Thus 
one may derive a Ps/n( z ) value for each point on the spec- 
trum once its S/N cont has been measured. At these faint 
limits, the signal to noise in the red becomes prohibitive 
when one encounters sky emission lines, and so we have cho- 
sen to apply this measurement to 7 windows between sky 
emission, at: 5000-5550A, 5600-5800A, 6000-6200A, 6600- 
6800A, 7100-7200A, 7600-7700A, & 8100-8270A. These cor- 
respond to redshift windows for [Oil] of z=0. 34-0. 49, 0.50- 
0.56,0.61-0.66,0.77-0.82,0.91-0.93,1.04-1.07 & 1.17-1.22. Fig- 
ure ^| shows 3 examples where Ps/n{ z ) has been calculated. 
To obtain a redshift constraint, we define a threshold for 
Ps/n{z). We can then use the point at which the interpo- 
lated Ps/n(z) curve rises above this threshold as a limit on 
the redshift. Here we choose a value of Ps/n(z) — 0.5 and 
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Figure 3. Redshift completeness of the sample versus mean sur- 
face brightness (/ifi, mea n), peak surface brightness (fiR^max), ob- 
served magnitude (R bserved) an d magnitude corrected for lens- 
ing amplification (R CO rrected)- Light, heavy and solid shading 
represent arclets for which redshifts were undetermined, tenta- 
tive and firm respectively. It is evident that our success rate for 
obtaining good redshifts falls as both the surface brightness and 
total observed magnitude fall. However, the bottom right panel 
indicates that, with the help of lensing amplification, we have 
obtained redshifts down to R ~ 25.5. 



quote the redshift limits obtained for these q = 3 spectra in 
Table |. 

There are several caveats to this process. Firstly, the 
value ps/n( z ) is only valid within the sky window where 
the signal to noise was measured. It is possible that a sig- 
nificant equivalent width of [Oil] could remain undetected 
between these windows if the continuum signal to noise in 
the sky band was low enough. Secondly, we assume here that 
the distribution of observed equivalent widths corresponds 
to those of the LDSS-2 deep redshift survey (Colless et al. 
1990). We note that this produces a conservative estimate 
since our survey pushes significantly deeper in magnitude 
and redshift thus making it likely that the true [Oil] equiv- 
alent width distribution of our sample contains more large 
equivalent widths than the one we have used. This would 
lead to an over estimate of Ps/n( z ) f° r an y given continuum 
signal to noise. Finally, we point out that this method will 
produce erroneous results if our sample does not correspond 
to the general field in its distribution of equivalent widths. 
To summarize, we have used our detection limit for emission 
lines to derive redshift limits for those spectra which show 
good continuum signal to noise but no features. 

Figure [?] shows how our redshift completeness varies 
with the target colours. It can be readily seen that more 
redshifts were obtained for blue objects than red. This prob- 
ably reflects the greater strength of [Oil] emission in the 
blue objects as hypothesised in our selection method. This 
is confirmed in Figure ^| which shows the correlation be- 



Figure 4. Continuum signal to noise ratio versus observed [Oil] 
equivalent width for the arclets with q = 1 (after Colless et al., 
1990). The curve shows a detection limit of S/Ni ine = 5. Due to 
the proximity of [Oil] to the A5577 sky emission line in #317, the 
equivalent width shown for this object is only a lower limit. Also 
shown is one candidate with q = 2 (#362) which falls just below 
the line. 



Figure 5. Examples of the 'p(z)' curves derived from the signal 
to noise arguments given in §2. The vertical arrows show how the 
redshift constraint is defined. The three examples show one case 
giving a single lower bound (#322), one giving both upper and 
lower limits (#492) and one giving no bound at all (#200). The 
horizontal bars at the bottom of the plot show the size of the 
redshift windows within which the probability value applies. 
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(B-I) 

Figure 6. The correlation between observed equivalent width 
and (B — I) colour for those targets where [Oil] was detected. The 
shaded area encloses the region occupied by the E/SO sequence 
of the cluster. [Oil] emission was detected in only one cluster 
member (#342) and in fact this object is well outside the cluster 
sequence. It is thus clear that the blue colour selection criterion 
proved a good predictor of the presence of [Oil] . 

tween equivalent width in [Oil] and (B — I) colour. Here, 
the shaded region marks the cluster sequence colour. [Oil] 
emission was detected for only one cluster galaxy which, in 
fact, lies well outside the cluster sequence. There is a clear 
anticorrelation between (B — I) colour and equivalent width 
such that the reddest galaxies have little or no [Oil] emis- 
sion. This confirms the results seen by Colless et al. (1990) 
and Glazebrook et al. (1995). We conclude that our colour 
selection was a useful way to increase the likelihood of de- 
tecting [Oil] emission (and thus obtaining a redshift), as well 
as a method for selecting non-cluster objects. 

Of the 37 arclet targets, we were able to obtain a total 
of 32 redshifts, 18 of which were of good quality (q = 1) 
and were behind the cluster (z s > 0.2). Figures ^[ and 
PH , show the extracted spectra, 2-dimensional spectra and 
zooms of the arclets, while the redshift catalogue is presented 
in Table |[ 

Several of the spectra merit individual discussion espe- 
cially with regard to their redshift determinations. First we 
deal with the multiple images and follow with those that are 
singly imaged. 

289 (R=20.5, z=1.034) & 359, (R=20.3, z=0.702) 
The spectra for these giant arcs are discussed by Pello et 
al. (1992). 

328 (11=22.0, z=0.702) & 389 (R=21.5, z=0.702) Both 
part of the #359 multiple system showing weak features 
corresponding to the #359 redshift of z — 0.702. #389 ap- 
pears seen through the disc of the spiral galaxy #373 and 
consequently the features are so weak that we assign q = 2. 

273 (R=21.9, z=0.800) The spectrum of this arc is heavily 
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(B-I) 

Figure 7. The (B — I) distribution of the spectroscopic sample. 
As before, light, heavy and solid shading represent undetermined, 
tentative and confirmed redshifts respectively. 

contaminated with light from the nearby elliptical #268 (see 
Figure ^). However, as may be seen from the 2-d spectrum, 
the emission line is very prominent at the lower edge of the 
#268 continuum. Only the 3 rows containing the emission 
line were extracted and no Gaussian weighting was used. 

384 (R=21.2, z=2.515) The spectrum of this spectacular 
arc has already been discussed in a recent letter (Ebbels et 
al. 1996.) Many rest frame ultra-violet metal lines are seen, 
along with Lyman-a in absorption at z — 2.515. 

444 (R=22.7, z=1.030) The spectrum has good contin- 
uum signal to noise and is extremely blue. We base our iden- 
tification on the Fell doublet at A2587 and A2600. This is in 
good agreement with the predicted redshift of z = 1.1 ±0.1. 
However this remains a tentative (q = 2) identification. 

145 (R=22.1, z=0.628) [Oil] is seen strongly in this spec- 
trum, along with some other more tentative features. 

158 (z=0.721) This spectrum could be a blend of 2 sources 
as may be inferred from Figure |^. We identify several fea- 
tures (principally [Oil]) at a redshift of z = 0.723 while other 
features, including the continuum, are more consistent with 
a redshift of z = 0.167. We suggest that most of the contin- 
uum comes from a bright cluster member while the fainter 
(by 0.7mag), diffuse object is a background galaxy (B=22.5, 
henceforth denoted #1581) and produces the [Oil] emission. 
Unfortunately the ground based seeing precludes any spa- 
tial identification of the emission line with one particular 
component. 

159 (R=21.9, z=0.564) [Oil] emission dominates this 
spectrum which also includes several absorption features. 

205 (R=20.3, z=0.693) This spectrum has good signal to 
noise with both absorption and emission features around the 
Balmer break visible. 

231 (R=22.0, z=0.563) The ground-based seeing causes 
this spectrum to be a blend of 231 and compact component. 
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These have R magnitudes of 22.0 and 23.4 respectively, so 
that no more than 21% of the light can be contributed by 
the compact component and it is unlikely that the emission 
line arises from it alone. 

238 (R=21.9, z=0.635) A single emission feature identi- 
fied with [Oil] dominates this spectrum. 

242 (R=21.8, z=654) Again, [Oil] emission is seen along 
with several absorption features at lower signal to noise. 

297 (R=22.4, z=0.450) [Oil] emission is very strong here 
and is accompanied by tentative identifications of the 
#/3/[Offl] triplet. 

306 (R=23.7, z=0.450) [Oil] emission along with several 
other more tentative identifications lead to the redshift for 
this arclet. 

317 (R=20.1, z=0.474) The [Oil] emission here falls very 
close to the A5577 sky line and so we can determine a lower 
limit only on its equivalent width. However, several other 
features confirm the redshift. 

381 (R=20.5, z=0.521) This spectrum, which has one of 
the best signal to noise ratios of all those obtained, exhibits 
many features from Mgll absorption to Hf3. 

431 (B=23.5, z=0.675) Despite its proximity to the A6300 
night sky emission line, the [Oil] emission in this arclet is 
strong and unmistakable on the 2-d spectrum. 

456 (R=22.2, z=0.538) Several absorption features in- 
cluding the K/H doublet are visible with good enough signal 
to noise to obtain a redshift for this arclet. 

464 (R=22.9, z=0.476) A strong emission feature identi- 
fied with [Oil] plus a possible G absorption band lead to a 
redshift for this arclet. 

467 (R=20.5, z=0.475) Excellent signal to noise provides 
a sure redshift for this arclet, with up to 14 identifications 
from [Oil] emission to A5268 absorption. 

We also mention the object #323 which, being very 
elongated but at an angle > 45° from the shear direction, 
would cause problems for the model were it at a redshift 
much higher than the cluster. Identification of several ab- 
sorption features in this spectrum reveal it to be a cluster 
member at a redshift of z — 0.179, thus posing no threat 
to the model. Finally, we note that of our new redshifts, 
7 arclets are found in three redshift structures. The sensi- 
tivity to small-scale redshift clustering is an artifact of the 
small angular extent of the survey and is a factor affecting 
all small-area surveys (including the photometric redshift 
surveys in the Hubble Deep Field, see §5). Thus, gaining 
a statistically-reliable view of the field redshift distribution 
will require inversion through a number of lenses. 

The 18 new redshifts bring the total number of arclets 
with good redshifts in Abell 2218 to 20 - making it the 
largest number for any lensing cluster to date. Table ^ lists 
the results in order quality with multiple image systems ap- 
pearing in the top 2 sections. The following analysis of the 
inversion method and our conclusions are based solely on the 
best quality (g = 1) arclets, though we expand the sample to 
include the q = 2 arclets where appropriate. For complete- 
ness we also list the q = 3 arclets and their redshift limits 
from the signal to noise argument above in order to give an 
indication of the sample which might be reached with a 10m 
class telescope. 



3 COMPARISON WITH LENSING 
PREDICTIONS 

We now compare the redshifts obtained with those predicted 
by the lens inversion method outlined in KESCS. 

The arcs whose redshifts can be predicted fall into two 
classes: multiple images and single images. The redshift pre- 
dictions for the multiple images have much smaller uncer- 
tainties than those for singly-imaged sources because of the 
larger number of constraints implied when several images 
of the same source are taken together. It is therefore par- 
ticularly important to prove that that the inversion works 
for these sources. However, in the case of a singly-imaged 
arclet, the uncertainty in the redshift prediction is more 
typical of that expected for the larger fainter population 
of arclets addressed by KESCS. The uncertainty here de- 
pends on the redshift, the size of the local shear and the 
unknown intrinsic shape of the source. The average uncer- 
tainties listed in KESCS were 8z/z = 22% for multiple im- 
ages and 8z/z = 30% for single images. 

In order to understand what level of agreement we ex- 
pect between our inverted redshifts and measured spectro- 
scopic ones, it is useful to consider the method by which the 
inverted redshifts are obtained. For any particular singly im- 
aged arclet at any point on our image, our mass model can 
predict magnitude and direction of the shear induced by the 
lens. The magnitude of this shear depends on the redshift 
of the arclet and increases from a minimum at the clus- 
ter redshift to an absolute maximum as the arclet redshift 
approaches infinity (This limit is actually replaced by a red- 
shift of ~ 4 in our model). The observed shape of the arclet 
is a combination of its intrinsic shape and this redshift de- 
pendent distortion. Conversely, knowing the observed shape, 
and given a redshift, one can predict the intrinsic shape us- 
ing the calculated shear. As we move the arclet in redshift 
from the distance of the cluster out to infinity, its predicted 
intrinsic shape will change, eventually tending to a constant 
shape at large redshift. 

We may represent the shape of any object on a plot of 
ellipticity vectors f as shown in Figure [12. Using the no- 
tation of KESCS, f = re 2ie and r = (a - fe 2 )/(2afe) with 
a and 6 the semi-major axes of the equivalent ellipse. The 
vectors may be plotted with respect to the axes of the CCD 
(large plots) or with respect to the axis of local shear for 
each arclet (insets). In these local shear axes, the lensing 
transformation conserves the y component of the elliptic- 
ity, T y and thus the process of moving the arclet in redshift 
simply corresponds to moving its point in the —x direction 
across the f-diagram. As we move the arclet out in redshift 
and thus to the left in t x across the plot, we can ask what is 
the likelihood that any given galaxy in the field would have 
the shape defined by this point on the f diagram. This like- 
lihood will correspond to the value of the ellipticity distri- 
bution of field galaxies, p(f) (represented by the contours on 
the plots) at that point. We can therefore define a redshift 
probability function p(z) for each arclet, corresponding to 
the likelihood that its intrinsic shape for any redshift would 
be found in the general population of field galaxies. The p(f) 
distribution peaks at an ellipticity of r = and has circular 
symmetry. Thus the maximum in our probability, p(r x ) or 
correspondingly, p(z) will occur when r x — 0. We therefore 
assign a most likely redshift z op t to each arclet defined by 
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Table 2. Arclct Spectroscopy in Abell 2218. We list first the multiple image systems, followed by the singly imaged ones. Only arclets 
with z s > 0.2 are listed, except in the case of #323 (see text). The arclets for which no redshift prediction was possible were added 
to the sample at points on the masks where suitable bright arclets with inverted redshifts were unavailable. COLUMN NOTES: (1) 
Rcorr calculated assuming z = z op t- (2) Quality (q): l=certain/probablc, 2=possible, 3=limit only (3) Prom WFPC-2 F702W image 
transformed to R as detailed in KESCS. Estimated error 0.2 mag. (4) For those objects with undetermined redshifts, we quote limits on 
the redshift derived from the argument in §2. COMMENTS: (1) Class III; No maximum in p(z): z op t — > oo (2) Class I; No maximum in 
p(z): z op t < 0.2 (3) Serendipitous object (4) Arclet off WF frame, photometry quoted in B band (5) Redshifts determined by Pello et 
al., 1992. 
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Figure 8. LDSS-2 spectra for arclets in Abell 2218, listed in order of identification number. All spectra with q = 1 are shown, while 
arclets #362, #389, & #730 are included for interest. (Arc #730 is shown with a rest frame wavelength scale corresponding to its 
predicted redshift.) See text for individual discussion. The spectra have been fluxed and smoothed to the effective resolution of LDSS-2, 
but are not corrected for atmospheric absorption. 



Figure 9. LDSS-2 spectra for arclets in Abell 2218. (Continued) 



the redshift which places it at t x — on the diagram of 
Figure 

This method divides the singly imaged arclet popula- 
tion into three classes, two of which preclude a determina- 
tion of z op t (see also Figure 7 of KESCS). The first, class 
I, corresponds to those arclets with t x < 0. These cannot 
be moved to r m = by placing them behind the cluster and 
must therefore be either cluster or foreground galaxies, or 
background galaxies which are not distorted enough to lie 
within a 45° cone about the axis of local shear. Therefore 
no meaningful z opt may be determined for them. The sec- 
ond, class II arclets constitute those which may be moved 
to t x = by placing them at a particular redshift behind 
the cluster and a determination of z opt is possible for these 
arclets. Finally, particularly in areas away from the cluster 
core where the shear is weak, we have some arclets (class III) 
which although aligned closely with the shear axis, do not 
reach t x — even at infinite redshift. Again, no meaningful 
z op t can be given for these arclets. 

We must emphasize that the redshift prediction z opt 
only has physical meaning in a statistical sense. That is, 
while the shape corresponding to z op t may not in fact be 
the actual intrinsic shape of the object, the distribution of 
these shapes across the population of inverted objects must 
be drawn from the distribution of field galaxy shapes, p(r). 
If this is true, then while the individual z opt values may not 
always correspond to the true redshift, the redshift proper- 
ties of the population as a whole (such as < z > and shape 
of the n(z)) should be correct. Later we assess the manner 
in which our inverted shape distribution (Figure |l^, left) 
samples the field galaxy p(f ) but first we discuss individual 
cases of inversion, starting with multiple images. 

In KESCS' model, seven multiply imaged systems were 
used to constrain the mass model. These were identified us- 
ing similarities in their HST morphologies, colours, surface 
brightnesses and on general assumptions about the lensing 
geometry. Only two of these (#359 and 289) already had 
redshifts from the spectroscopic study of Pello et al. (1992). 
Furthermore, none had actually been confirmed spectroscop- 
ically as multiple. Thus our redshift for #328, a compo- 
nent of the 5-image #359 system is a strong vindication 
of the mass model and the ability of HST to identify mul- 
tiple images. The spectrum of #389 is consistent with its 
being a third image of is system although highly contami- 
nated by light from the obscuring cluster spiral #373. In- 
deed, the identification of this as a 5-image system was the 
major motivation of KESCS to extend their model to in- 
clude galaxy scale mass components. These 5 images could 
not be modeled using cluster scale mass components alone. 
Therefore, this confirmation reinforces our confidence in the 
mass model down to very small scales (~ 75fcpc) . The #384 
system consists of one merging image pair (arc #384) and 
its counter- image (arc #468). Ebbels et al. (1996) confirmed 
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Figure 13. The fractional difference between zi en3 i ng and z sp ec 
for the A2218 sample. Single images are represented by squares 
and multiple images by triangles. The open triangles show mul- 
tiple images that were used as constraints on the model. Objects 
with q = 1 are shown with solid error bars, while those with q = 2 
have dashed error bars. 



that both halves of the merging pair have the same redshift 
(z — 2.515) and also found the spectrum of #468 to be 
compatible with this redshift, although its signal to noise is 
poor. Of particular significance, and also discussed briefly in 
Ebbels et al., is the fact that the spectroscopic redshift of 
this system is in excellent agreement with the lensing pre- 
diction: z = 2.8±q I (KESCS). The redshift of this system 
is particularly important due to its proximity to a critical 
line, meaning that the predicted redshift is extremely sensi- 
tive to changes in the mass model. In summary, therefore, 
four images from two multiple systems are consistent with 
the earlier mass model and the most highly-magnified has a 
redshift within the expected uncertainties. 

Turning to the singly-imaged arclets, Figure |l^ shows 
the fractional difference between z op t and z spec for the arcs 
and arclets where both a spectroscopic and inversion redshift 
above z = 0.2 are available. Here the error bars refer to 
the 80% likelihood values discussed by KESCS (henceforth 
denoted by <7 80 %) and we include arclets with q = 2 for 
purposes of illustration. Although the sample of q — 1 arclets 
is small, it is useful to ask at what level of significance do 
the results agree with the predictions. 
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Figure 10. LDSS-2 spectra for arclets in Abell 2218. (Continued) 
Figure 11. LDSS-2 spectra for arclets in Abell 2218. (Continued) 



The arclets can be grouped according to the degree 
of success with which we predicted their redshifts. In total 
there are 11 arclets of class II and 9 = 1 for which mean- 
ingful redshift predictions can be made. Of these, 9 match 
their spectroscopic redshifts within 2a so %. As tests of the 
method, the best of these are arclets #231, #297, #317, 
#381 #467 & #456 which have small error bars leaving lit- 
tle likelihood that the predictive success could be a product 
of chance alone. The remaining 3 arclets (#145, #1581 & 
#238) have larger error bars due to their relatively large 
values of \r y \, but still agree tolerably well with their spec- 
troscopic redshifts. Finally there are two arclets whose pre- 
dictions do not match their spectroscopic redshifts within 
2cr g0% - objects #205 and #464. Morphologically, #205 is 
clearly a spiral galaxy (Figure [|) and this apparent lack of 
distortion so close to the cluster core leads to the low redshift 
prediction. #464, however appears quite distorted in the di- 
rection of the local shear and its predicted redshift is corre- 
spondingly high. The disagreement between prediction and 
measurement for the arclets #205 and #464 must therefore 
result from the natural dispersion of image shapes on the 
sky, leading to chance alignments in the source plane away 
from and along the shear direction respectively. Given the 
shape of the field galaxy ellipticity distribution as illustrated 
by contours in Figure jl2|, it is not unusual to find galaxies in 
these regions (reasonably large \t x \ but small \t v \). Thus we 
conclude that these two cases represent merely a statistical 
fluctuation in this small sample. 

We can now consider, in the light of Figure [l3|, the na- 
ture of the error in the predictions for the population as a 
whole. Taking only the q = 1 arclets, the distribution shown 
in Figure [l3| has a mean of < 5 z /z >= 0.07 and dispersion 
°"<Sz/z = 0.48, while including the q = 2 arclets brings these 
values to < S z /z >= 0.09 and <t<s,/ z = 0.41. We therefore 
expect the fractional error in the predicted quantities of the 
whole population (such as < z >) to be of order 10% or less, 
despite the dispersion of individual points. 

The conclusion from the above discussion is that while 
z opt and z spec agree well for many of the arclets, we cannot 
expect this to be the case for all arclets. This is simply due to 
the unknown position of f within the natural distribution 
of source shapes. As mentioned above, it is useful to con- 
sider how the predicted ellipticity distribution of our sample 
compares with that seen in the field generally. Returning to 
Figure the lower two panels illustrate how both the spec- 
troscopic and fainter samples trace the field galaxy distribu- 
tion (contours, measured from ~ 10000 galaxies between 
I = 18 - 25 taken from HST WFPC2 fields in the MDS and 
Groth Strip Surveys (Ebbels et al. in prep.)). Two effects can 
clearly be seen. The first is the tendency in the image plane, 
for the spectroscopic sample (filled symbols) to be more ex- 
tended than the contours. This justifies our selection of these 
objects as arclet candidates - they are more extended than 
the field population in general. The second effect is a re- 
sult of the method itself, specifically that our choice of the 
optimum redshift of an object will correspond to its low- 



est obtainable ellipticity. Thus the distribution of points in 
the source plane will tend to be more centrally concentrated 
than in the field and this can readily be seen in the Figure. 
Notwithstanding these biases however, it is clear that our 
predicted distribution favors no particular orientation and 
is centrally peaked - the two features most noticeable about 
the unlensed field distribution. 

In using this technique to probe the redshift distribu- 
tion of faint galaxies, we are mainly concerned with certain 
statistical properties of this distribution such as the mean 
< z > and the dispersion a z . We are therefore compelled 
to compare the predicted values of these quantities with 
their actual values for our sample. We show the redshift 
distributions derived from the inversion method and from 
spectroscopy in Figure [w], where we have ensured that both 
distributions correspond to that subset of our sample pos- 
sessing both inversion and spectroscopic redshifts. Clearly 
the two distributions are similar in form with the bulk of 
the arclets below a redshift of one. Using just those arclets 
with q — 1, we find < z 3p( , c >= 0.68 while the predicted value 
is < z op t >= 0.69, a difference of just 1.5%. Expanding the 
sample to include the q — 2 arclets, we find < z spec >= 0.63 
and < Zopt >= 0.61, a difference of 2%. In this case therefore, 
the mean redshift is very accurately predicted by the tech- 
nique. Looking at the shape of the distributions, one may 
compare their dispersions. Considering the q = 1 data, we 
find cr z c = 0.49 and cr Zopt = 0.62 while including the data 
with q = 2 gives <r s „ ae = 0.42 and u Zopt — 0.53. Bootstrap 
resampling of the true (spectroscopic) distribution indicates 
that the widths of the inversion and spectroscopic distri- 
butions are not statistically distinguishable at greater than 
~ 75% confidence. We do expect, however, that the uncer- 
tainties in the inversion redshifts will increase the predicted 
dispersion over the spectroscopic one. We may therefore con- 
sider the inversion-predicted dispersion better as an upper 
limit on the true width of the redshift distribution. In sum- 
mary then, this analysis gives us confidence that our method 
can find the correct mean redshift for an unknown distribu- 
tion, but that the width of the distribution is less well de- 
termined. For the mean redshift, the agreement between the 
predicted and spectroscopic values is especially encouraging 
when one considers the small sample size. 



4 IMPLICATIONS FOR THE MASS MODEL 

In this section we return to the mass model used by KESCS 
and investigate to what extent the new spectroscopic results 
further constrain it. 

First, we consider only the spectroscopic redshifts of the 
multiple images, namely: #359 (z=0.702), #259 (z=1.034), 
#384 (z=2.515), and #444 (z= 1.030) and optimize the dif- 
ferent cluster and galaxy scale components following the 
KESCS prescription. Following the analysis of Natarajan 
& Kneib (1997), we also vary the exponent a of the scal- 
ing law of the truncation radius (r cu t = r* ut (L/L*) x ' cl ). We 
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Figure 12. The lens inversion method. We show plots of the f ellipticity vectors for the arclets in the source plane (left panel) and 
the image plane (right panel). For each plane, we show the f vectors with respect to the frame of the local shear (upper panel) and the 
CCD frame (lower panel). In all plots we show linearly spaced contours of the field galaxy f distribution as measured from HST field 
survey images (Ebbels et al. in preparation). Solid symbols represent the spectroscopic sample while open ones indicate the full sample of 
objects with areas over 50 pixels (137 arclets.) The inversion method is clearly seen at work in the upper panel. In the image plane (top 
right), most of the arclets have t x > since most of them (especially the spectroscopic sample) arc aligned with the local shear. Placing 
each source at a redshift behind the cluster translates the corresponding point on this diagram in x only. For each arclet we choose that 
redshift at which r is minimized, and thus the invertible arclets (class II) will end up on the line t x = in the source plane (top left). 
In the source plane, class I arclets (not sufficiently distorted to invert) lie in the region t x < and class III arclets (too elongated to 
invert) lie in the region t x > as is clearly seen. Moving to the CCD frame, one can see that the distribution of points in the image 
plane (bottom right) is more extended than that of field galaxies in general (justifying their selection as elongated arclets.) In the source 
plane (bottom left) the effect is reversed since we choose the roundest shape for each object and the distribution of points is more peaked 
than the field. It is clear that when referred to the CCD coordinates, there is no preferential direction, even in the source plane. 
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Figure 14. The distribution of our sample with spectroscopic 
(left) and inversion (right) redshifts. Objects with firm (q = 1) and 
tentative (q = 2) redshifts are denoted by solid and heavy shading 
shading respectively. Here we have ensured that both histograms 
correspond to the same sample of objects which include multiple 
images (except those used as model constraints). In each plot 
the point with error bar marks the mean and dispersion of the 
distribution for the q = 1 sample. 



tried two different values: a = 4 to mimic a constant mass 
to light ratio for the galaxies, and a = 2.5 to mimic a 'Ko- 
rmendy' relation between the effective radius of the mass 
and the luminosity. The resulting best fit mass distribution 
does not differ greatly from the one presented in KESCS. 
Less mass is required between the two main clumps than in 
the KESCS model, and the a = 2.5 exponent is preferred 
leading to a mass-to-light ratio of 9 for an L» galaxy. 

The improved model in turn allows us to put better 
constraints on the redshifts of the other multiple images. 
In particular, the #Hl-3 multiple arc appears to be clearly 
different from the arclet #273 (z=0.80) as its redshift is now 
predicted to be z — 1.9 ± 0.2. In fact, closer inspection of 
the HST/WFPC2 image reveals a small offset between the 
#273 and #Hl-3 systems of ~1 pixel, suggesting that these 
images are indeed separate systems. The multiple system 
#730 suffers the same translation to higher z as HI, 2, 3 
and is now expected to be at a higher redshift than ~3. 
(With the relatively low redshift of A2218, our ability to 
state tighter constraints for high redshift arcs is hindered by 
the weak dependance on redshift of the Dls/Ds ratio.) 

We then considered both the singly and multiply im- 
aged arclets whose spectroscopic redshifts are presented 
here. Their images were mapped back to the source plane 
to compare their intrinsic shape distribution to that of the 
unlensed field (shown by contours in Figure [li]). These in- 
trinsic ellipticities form a circularly symmetric distribution 
in the r^-plane and, with few exceptions, have ellipticities 
within t < 0.2. This area contains over half of the entire un- 



lensed distribution and therefore the predicted shapes of our 
arclets well match those expected for unlensed field galaxies. 
To include the single image redshifts in the optimisation, we 
reran the optimisation procedure, adding to the x 2 sum a 
term of the form (rs/a T s) 2 for each arclet, where a T s is the 
width of the intrinsic ellipticity distribution. The difference 
between the two optimizations is minor, demonstrating that 
we are now very close to the true mass distribution within 
the boundary of the HST/WFPC2 field. 

Finally, we used the updated model on the catalogue of 
faint arclets with areas above 50 pixels (as used in KESCS) 
to compute the mean redshift versus magnitude prediction. 
This is shown by the thick line on Figures |l^ and [l(] where 
we compare it to the prediction of KESCS (thin lines). 
Clearly, while there may be some differences in z op t for indi- 
vidual arclets, the difference in the prediction for the whole 
sample (even in the faintest bin) is small ( <5%). The new 
mass model is the best we can achieve with the present red- 
shift constraints. 

We now discuss how the residual uncertainties in the 
mass model can affect our redshift predictions. Of course, 
we can only change the mass model in such a way that it 
still reproduces the multiple imaging seen in the HST im- 
age. In the modelling process (see KESCS) we use a \ 2 sum 
to represent how well the predicted parameters of the mul- 
tiple images match those observed. Thus, in perturbing our 
model, we will try to keep the resulting x 2 as close to that 
of the fiducial model as possible. A number of tests may be 
made, corresponding to uncertainties in the model on both 
small and large scales. 

Firstly, we may investigate the effect of the galaxy scale 
components. These were included in the model solely to 
reproduce several of the multiple image systems correctly. 
Only a few were actually needed to do this but a larger 
number (~ 30) were added corresponding to a given mag- 
nitude limit in the interests of consistency. We tested their 
effect on the redshift predictions by removing them at ran- 
dom from the model and predicting the redshift - magnitude 
relations, as shown in Figure The triangles, circles and 
squares show the mean and dispersion in < z > for each bin 
when clumps were removed in groups of 1, 2 and 5 respec- 
tively. We also show the predictions of our new model (thick 
line) and that published by KESCS (thin line) where both 
include the statistical correction for cluster and foreground 
contamination as described in KESCS. The only individ- 
ual cases deviating significantly from the fiducial predictions 
were those where the galaxy scale clump corresponding to 
the cD galaxy was removed. However, since an external con- 
straint in the form of the cD velocity dispersion exists for 
this clump, we do not consider its removal to be a valid 
test and exclude these cases from the analysis. In all three 
cases, the dispersion in < z > for each bin never rises above 
cr< z > = 0.09 and the average fractional dispersion is only 
<7< z > / < z >= 5.6%. Clearly, the mean redshift and its trend 
with magnitude are insensitive to the details of the mass 
model as represented by the galaxy scale clumps. 

On large scales (~ IMpc), the uncertainties in the mass 
model derive mainly from the small angular size of the HST 
image. Outside the region of multiple imaging, we obtain less 
information on the mass profile, and features such as the ex- 
tent of the main cluster clump or the presence of external 
mass clumps are poorly constrained. However the effect of 
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Figure 15. The effect on the redshift predictions of removing 
galaxy scale clumps from the model . Points with error bars show 
the mean and dispersion in < z > for each bin. triangles: One 
clump removed at a time, squares: 50 random pairs of clumps re- 
moved one at a time, circles: 50 random sets of 5 clumps removed 
one at a time. The three sets of points are offset in R for clarity. 
In each case, the thin line shows the mean redshift-magnitudc 
relation for the fiducial model (KESCS) while the thick solid line 
shows the relation predicted by our new model developed in §4. 



these uncertainties can be readily tested in a similar manner 
to that of the galaxy scale clumps. Figure ^ shows the effect 
on the magnitude-redshift relation of two such tests: adding 
external large scale clumps outside the HST area (squares) 
and varying the cut radius of the main cluster clump (cir- 
cles). Again, each point shows the mean and dispersion in 
< z > for each bin, and we show the predictions of KESCS 
and our new model by thin and thick lines respectively. For 
the first test, large scale mass clumps varying in velocity dis- 
persion (500 — 2000fcm/s) and core radii (700 - 1500fcpcf] 
were added to the model at the same redshift as the cluster 
and four different positions 1.9 Mpc from the centre of the 
cluster (outside the HST frame). 48 different combinations of 
clump parameters were tried and in each case the model was 
re-optimized in order to reproduce the multiple image con- 
straints as well as possible. Only the models giving a residual 
X 2 within 30% of the fiducial model value were considered. 
Again, the maximum dispersion is low, a <7> = 0.10 and the 
average fractional dispersion is just <r< z > / <z >= 6.4%. For 
the second test, the cut radius of the main cluster clump, 
r cu t, was varied between 400 — 2000fcpc, re-optimizing the 
model to reproduce the multiple images as before. Here the 
maximum dispersion is <r< z > = 0.04 while the average frac- 
tional dispersion is a <z> / <z> — 4.5%. Once more, we find 



t Here and in the following we assume H = 50 km s 1 Mpc 1 
and q = 0.5 



Figure 16. The effect of adding large scale external mass concen- 
trations to the model {squares) and of varying the cut radius of 
the main cluster clump [circles). Each point and error bar shows 
the the mean and dispersion in the < z > for each bin and the 
points are separated in R for clarity as before. Again the thin 
solid line represents the fiducial model of KESCS while the thick 
line shows the predictions from our new model of §4. 



that the uncertain features of our model on large scales do 
not affect the results appreciably. In summary then, our red- 
shift predictions for the faint population behind A2218 are 
unlikely to be much affected by the remaining uncertainties 
in our mass model. 



5 PROPERTIES OF THE FAINT POPULATION 

The results of our spectroscopic program presented in sec- 
tion §3 give strong support to the inverted redshift distribu- 
tion derived for fainter arclets viewed through Abell 2218. 
Moreover, the discussion of §4 demonstrates that further 
changes to the mass model permitted by our spectroscopic 
data will not greatly affect these redshift predictions. 

For arclets with well measured shapes, such as our spec- 
troscopic sample, we found in §3 an error in the mean red- 
shift of ~ 2%. Using fainter arclets with more poorly mea- 
sured shapes would increase the errors for each individual 
arclet and thus increase the uncertainty in the mean red- 
shift for the population as a whole. This point was addressed 
fully in KESCS where it was shown using simulations that 
reliable shapes could be derived for arclets with isophotal 
areas above 50 pixels. We must therefore restrict ourselves 
to this area limit in any analysis of the population. The re- 
sults (Figures 15-16) show a mean redshift < z >~0. 8-1.0 
for amenable sources selected with R ~25.5 (B ~ 26-27). 

A low mean redshift at such a faint limit has consis- 
tently emerged from statistical lensing studies (Smail et al. 
1994, KESCS). Recently, Luppino & Kaiser (1997) have re- 
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ported the detection of a significant weak shear signal from 
the bluest half of the faint field population behind a distant 
cluster, inferring a high mean redshift (z >> 1) for this sub- 
sample. However, we note that by combining the blue and 
red subsets from Luppino & Kaiser's analysis the mean red- 
shift they would infer for the whole faint galaxy population 
would be much closer to the no evolution expectation (in line 
with the results presented here) especially when allowance 
is made for the various uncertainties in their measurements 
and analysis. 

The only independent technique which has been used to 
investigate the redshift distribution of galaxies beyond the 
spectroscopic limit is that based on multicolour photometry 
(the photometric redshift technique). The popularity of this 
technique has increased considerably of late because of pre- 
cise data available for the Hubble Deep Field (HDF). It is 
thus of interest to compare our results with those of work- 
ers who have determined photometric redshifts in the HDF. 
Several surveys have been published using the HDF opti- 
cal photometry combined with spectral templates derived 
either from population synthesis or local spectral energy dis- 
tributions or both. The advantages of adding near infrared 
photometry to follow the usual spectral features to redshifts 
beyond ~ 1.3 are clear, but are difficult to compare with our 
results since they are usually magnitude limited in some in- 
frared band. Here we compare our redshift distribution with 
the results of Mobasher et al. (1996) and Lanzetta et al. 
(1996). 

Such comparisons are, however, hindered by the differ- 
ent techniques used and, in particular, that it is difficult to 
construct a magnitude limited sample from the lensing in- 
version method. Firstly, since the magnification of images 
varies across the cluster, our observed magnitude limit does 
not correspond to a sharp intrinsic magnitude limit, rather 
the lensing sample suffers from a more gradual fall off in 
counts with magnitude. One approach would therefore be 
to choose an effective magnitude limit corresponding to the 
point where the counts fall significantly below those in the 
field. The second complication arises from the area limit dis- 
cussed above. In this case, we can choose to select the same 
proportion of faint galaxies as a function of area from both 
the lensing and photometric redshift samples. These choices 
correspond to first making a magnitude cut at R ~ 24 
(/ ~ 23.5) and then taking only the largest 59% of the sam- 
ple at this limit. Figure |l7| shows the comparisons, where 
we have chosen to use a slightly deeper magnitude limit 
(/ < 24) for the photometric redshift distributions in order 
to obtain a comparable number of galaxies from the HDF. 

It can readily be seen that the agreement is much bet- 
ter when we compare to the results of Mobasher et al. than 
Lanzetta et al. (although the numbers of galaxies are small.) 
Excluding the galaxies below z — 0.2 (which cannot be 
probed by our lensing analysis), we find the lensing n(z) 
has a median of z — 0.4 with only 7% of objects appear- 
ing above z = 1.5 while Mobasher et al.'s distribution has a 
slightly higher median redshift (z = 0.5) but with a longer 
tail (19% above z = 1.5). The corresponding numbers for 
the Lanzetta et al. distribution are z = 0.8 and 18%. How- 
ever, as one goes fainter, the agreement worsens. Figure |ri| 
shows the predicted redshift distributions for our faintest 
magnitude slice 24 < R < 25.5. At this limit, the lensing 
method predicts a low median redshift of z = 0.7 while the 
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Figure 17. Comparison of the lensing n(z) (shaded) with 
the photometric n(z) from the Hubble Deep Field analyses of 
Mobasher et al (top) and Lanzetta et al. (bottom) for the mag- 
nitude limit R < 24. Points with errorbars show the median of 
each distribution, with errors calculated from bootstrap resam- 
pling. In order to make a valid comparison, only the largest 59% 
of objects were used. It is clear that there is much closer agree- 
ment with Mobasher et al than Lanzetta et al, emphasizing the 
discrepancies between the two photometric n(z) determinations. 



corresponding values from optical multi-colour photometry 
are z = 0.9 (Mobasher et al.) and z — 1.1 (Lanzetta et al.) 
respectively. 

At the faintest limits, it is important to bear in mind 
that the sensitivity of photometric redshifts in the 1 < z < 2 
region is expected to be very poor as there are no significant 
continuum features at optical wavelengths. For z >1.3, the 
4000 A break disappears from the optical window and the 
Lyman break does not enter until z > 2. Is it therefore 
possible that some galaxies allocated to 1.5 < z < 2.0 by the 
photometric technique may actually lie at lower redshifts ? 

An important recent development in this regard is the 
addition of JHK ground-based near-infrared data to comple- 
ment the earlier optical photometric redshift studies. Con- 
nolly et al. (in prep.) report a redshift distribution based 
on 4-colour optical plus 3-colour near infrared data for a 
J <23.5 sample which has < z >=0.9±0.1. This corresponds 
to a limit somewhat deeper than R = 25.5 but is a good 
indication of the changes that we might expect when near- 
infrared data is added to the studies discussed above. 



6 CONCLUSIONS 

We summarize our conclusions as follows: 

• We have conducted one of the deepest spectroscopic 
surveys attempted on a 4m class telescope and certainly the 
faintest with LDSS-2. The faintest redshifts were obtained 
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Figure 18. Comparison of the lensing n(z) (shaded) with the 
photometric n(z) of Mobasher et al. (top) and Lanzetta et al 
(bottom) for the magnitude slice 24 < R < 25.5 (23.5 < I < 25). 
Again, points with errorbars show the median of each distribu- 
tion, and only the largest 59% of objects used. The lens inversion 
method predicts a significantly lower median redshift than those 
based on optical multi-colour photometry at these limits. 



for sources with R—24 whose unlensed magnitudes reach 
R—25 in some cases. The major limiting factor is the effec- 
tiveness of sky subtraction and our limit is more appropri- 
ately constrained by surface brightness than integrated mag- 
nitude. We have achieved a useful success rate to hb. = 24 
arcsec -2 with completeness falling below 50% at hr — 23 
arcsec -2 . 

• The results of our spectroscopic survey show that the 
lensing inversion technique is reliable. The mean redshift 
of our spectroscopic sample determined from both inver- 
sion and spectroscopic redshifts agree to within ~ 2%. Since 
the inversion method is purely geometric in operation, this 
agreement gives us confidence when estimating the mean 
redshift of the larger population of faint arclets provided 
their individual shapes can be adequately measured. 

• We have included the new spectroscopic redshifts to 
further constrain our mass model and find that the redshift 
predictions change by a maximum of ~ 5%. Moreover, the 
predicted redshifts for a population of arclets with R ~25.5 
are insensitive to further changes in the mass model allowed 
by the multiple image constraints. We conclude that errors 
in the inferred mean redshift deriving from uncertainties in 
the mass distribution can remain only below the level of 

&<£>/ <2> <10%. 

• The new predictions for the faint population do not dif- 
fer significantly from those presented by KESCS. The mean 
redshift for the population of arclets amenable to inversion 
with R ~25.5 is broadly consistent with the no evolution 
prediction; we find <z>~ 0.8 — 1.0. 

• To a limit set by our magnitude and area cuts, the in- 



version prediction agrees reasonably well with results pre- 
sented by other authors using photometric redshifts deter- 
mined from data on the Hubble Deep Field. A direct com- 
parison is difficult since the inversion technique requires an 
area limited survey in contrast to one that is magnitude lim- 
ited. The agreement may be much better when near-infrared 
photometry is used in addition to optical colours. This can 
be understood via a greater sensitivity to the 4000 A break 
which is redshifted into the near-infrared at redshifts of in- 
terest. 

Application of this technique to other clusters with sim- 
ilar data has already begun (Kneib et al. in prep.). This will 
generate a much larger sample of inverted arclets and hence 
beat down the statistical uncertainties, as well as overcoming 
the effects of redshift clustering in the background popula- 
tion. Using cluster lenses over a wide range of redshifts will 
reduce the intrinsic uncertainties coming from the redshift 
dependence of the distance ratio, while extending the tech- 
nique into the near-infrared using NICMOS will allow us to 
select 1 < z < 3 galaxies in their restframe optical, to pro- 
vide a more robust measure of the star-formation density 
at these epochs. By expanding our data set in this way, wc 
shall be able to use this proven technique to derive robust 
and useful parameters for the faint population of galaxies at 
depths presently unobtainable by conventional means. 
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